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ABSTRACT 

We study the force balance and resulting acceleration of gas in general relativity basing on simulations of 
accretion on a stellar-mass, non-rotating black hole. We compare properties of acceleration in an optically thin, 
radiatively inefficient disk, and in an optically thick, super-critical disk accreting at 10 times the Eddington 
rate. We study both the average forces acting at given location and forces acting on a gas parcel along its 
trajectory. We show that the acceleration is not a continuous process - in most cases gas is accelerated only 
in short-lasting episodes. We find that in the case of optically thin disks gas is pushed out by magnetic field 
in the polar region and by thermal pressure and centrifugal force below the disk surface. In case of optically 
thick, radiative accretion, it is the radiation pressure which accelerates the gas in the polar funnel and which 
compensates (together with the centrifugal force) the gravity in the bulk of the disk. We also show that the 
Newtonian formulae for the forces are inadequate in the innermost and in the highly magnetized regions. 

Keywords: accretion, accretion discs - black hole physics - relativistic processes - methods: numerical - 
galaxies: jets 

outflowing gas along the line of sight have been observed 
for years both in case of accreting stellar mass compact ob¬ 
jects and active galactic nuclei (AGN) (King et al. 2013). 
It has been shown that winds are most prevalent in the soft 
state of black hole (BH) X-ray binaries (Ponti et al. 2012). 
Their presence is anticorrelated with the presence of radio 
jets (Miller et al. 2008; Neilsen & Lee 2009). They are highly 
variable (e.g., Neilsen et al. 2011; Neilsen & Homan 2012; 
Neilsen 2013) and carry matter out at significant rates (e.g., 
Ueda et al. 2004; Kallman et al. 2009; Miller et al. 2016). 

Deriving the properties of the outflow from the observed 
absorption lines is not straightforward. The observed ab¬ 
sorption occurs along the line of sight as a result of in¬ 
terplay between gas density and ionization, often far from 
the region where the outfiow was generated. In addition, 
the absorption profiles are superimposed on the relativisti- 
cally broadened reemission features coming from larger vol¬ 
ume (Miller et al. 2015). There is no consensus on whether 
the observed outflows in X-ray binaries are driven magneto- 
centrifugally (Miller et al. 2006, 2008; Kallman et al. 2009; 
Neilsen et al. 2016; Miller et al. 2016), thermally (Ueda et 
al. 2009; Neilsen et al. 2011, 2014; Diaz Trigo et al. 2014), 
or whether both mechanisms work at the same time (Neilsen 


1. INTRODUCTION 

Gas accreting on compact objects liberates large amounts 
of its binding energy. Even small fraction of this energy con¬ 
verted into radiation makes accreting systems extremely lu¬ 
minous. But radiation is not the only product of accretion. 
There is a growing evidence that accretion inevitably gen¬ 
erates outflows carrying outward significant amount of gas - 
sometimes much more than what reaches the compact object. 

Gas may be driven out of the accretion fiows in a couple 
of ways. Large scale magnetic fields may lead to magneto¬ 
centrifugal acceleration of gas along rapidly rotating field 
lines (Blandford & Payne 1982; S^dowski & Sikora 2010). 
Radiation generated by the infalling gas may exert outward 
pressure that accelerates the gas (Shakura & Sunyaev 1973). 
Einally, illumination of the outer regions by hot radiation 
coming from the inner region may drive thermal outfiow 
(Begelman et al. 1983). 

Blueshifted absorption features indicating the presence of 
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& Homan 2012). Active galactic nuclei show significant out- 
fiows as well (King et al. 201 1, 2012), some of which can be 
mildly relativistic (King & Pounds 2015; Chartas et al. 2002; 
Tombesi et al. 2010), providing physical mechanism for con¬ 
necting the SMBH and the properties of its host. 

Analytical modeling of accretion flows usually takes 
strong assumptions, e.g., by neglecting the outflows com¬ 
pletely (Narayan & Yi 1994) or prescribing their radial de¬ 
pendence (Blandford & Begelman 1999), and does not al¬ 
low to study their generation consistently. Multi-dimensional 
simulations are necessary for this purpose. However, they are 
limited, and a single simulation is not able to cover the whole 
accretion fiow. Woods et al. (1996), Luketic et al. (2010) and 
Higginbottom & Proga (2015) focused on the outer regions 
and studied the generation of thermal winds from the pho¬ 
tosphere irradiated by strong X-rays produced in the inner¬ 
most part. Other groups (e.g., Proga et al. 2000; Gammie 
et al. 2003; Ohsuga et al. 2009; Sadowski et al. 2014) stud¬ 
ied the innermost region where magnetocentrifugal and ra¬ 
diative driving operates and showed that, indeed, singificant 
outfiows emerge. 

In this paper we study the acceleration mechanisms acting 
on outflowing gas in optically thin and optically thick accre¬ 
tion flows by analyzing simulations of the innermost region 
performed with general relativistic radiation magnetohydro- 
dynamical (GRRMHD) code KORAL. To model these phe¬ 
nomena, we break down the individual forces acting upon 
the gas and determine the main acceleration mechanisms in 
each case. 

Similar modelling has been done previously by Yuan 
et al. (2015) and Takahashi & Ohsuga (2015), but this pa¬ 
per presents the first general relativistic modeling and di¬ 
rectly compares driving mechanisms in optically thin and 
thick disks performed within the same computational frame¬ 
work. We study only accretion fiows on non-rotating, stellar 
mass BHs, and limit ourselves to accretion with weak, not- 
saturated magnetic fields (SANE, in the formalism put for¬ 
ward in Narayan et al. 2012). 

The paper is organized as follows. In Section 2, we decom¬ 
pose the forces acting on a particle. Section 3 describes the 
simulation environment and the initial conditions. Sections 4 
presents the results and Section 5 discusses them and gives 
concluding remarks. 


1.1. Definitions 

We adopt the following definition of the Eddington accre¬ 
tion rate. 




-^Edd 


( 1 ) 


where 


% = 1 - a/1 


3^isco ’ 


( 2 ) 


is the efficiency of a thin disk (Shakura & Sunyaev 1973; 
Novikov & Thorne 1973) extending down to the marginally 
stable orbit located at Risco- According to this definition, a 
thin disk accreting at MEdd emits the Eddington luminosity. 
Eor zero BH spin, MEdd = 2.48 x IO^^M/Mq g/s. 

We use the term outfiow to denote any consistent motion 
with positive radial velocity This definition includes both 
the real outfiow and the turbulent outfiow. While the former 
eventually reaches infinity, the latter eventually rejoins the 
accretion fiow due to the turbulent gas motion in the disk To 
determine whether or not a fiow is part of the real outfiow, 
we define the relativistic Bernoulli parameter Be (Sadowski 
& Narayan 2015), 


Be = - 


T] pu^ 

pu^ 


( 3 ) 


where T,R,p and u^ are gas and radiation stress-energy ten¬ 
sors, gas density and four-velocity, respectively. The real out¬ 
fiow will satisfy Be > 0. However, it is important to note that 
the Bernoulli parameter of a fiow will often fluctuate. Despite 
having a negative parameter when it reaches the simulation 
boundaries, the fiow might still be part of the real outfiow, 
if it ultimately reaches Be > 0. This can typically be seen 
as a steady increase in the parameter over the course of the 
trajectory (Yuan et al. 2015). 

We define polar funnel as the conical region of low density 
between the disk surface and the polar axis. 

Hereafter, we use the gravitational radius rg = GM^ulc^ 
as the unit of length, and rg/c as the unit of time, and often 
assume c = 1. We also always assume that the underlying 
spacetime metric is fixed, i.e., is not perturbed by the accre¬ 
tion fiow. 


2 . EORCE DECOMPOSITION 
2.1. Conservation laws 

The three general conservation laws, that is conservation 
of mass, energy, and angular momentum, are described by 
the following set of equations. 


= 0, 

(4) 

in),, = Gy, 

(5) 

> 

1 

II 

s 

(6) 


where p is the gas density in the comoving frame, u^ = 
(u\ u^^ is the angular four-velocity, semicolon indicates the 
covariant derivative, and Ty is the magnetohydrodynamic 
stress-energy tensor, 

Ty = wu^Uy -r ]^b^6y - If by. (7) 

w = p Ug pg b^ denotes relativistic enthalpy, Ug and 
p^ = (y - i)ug are the internal energy and gas pressure, re¬ 
spectively, 7 = 5/3 is the adiabatic index. If is the magnetic 
field four-vector and Sy is the Kronecker delta. 
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Ry is the stress-energy tensor of radiation and Gy is the 
radiation four-force describing the interaction between gas 
and radiation defined by, 

( 8 ) 

where Gq and refiect the energy and momentum trans¬ 

fer due to absorption and scattering, and Comptonization, re¬ 
spectively. The former is given by (Sadowski et al. 2015), 

Gq = -p{Ka + Ks)R^^Uy - p (KgR^^UaUp + Ka^TlB^ , (9) 


yielding, 

di( ^Pgwu‘uy) = - dii ^pgp^6‘y) - di I 

+ dii ^b'by) + + V=gGv. 

(17) 

The left hand side of Equation (17) equals, 

di{ ^f^WU^Uy) = Uy • di{ • di(Uy) 

= Uy • di( ^[^{w - p)u^) + • diiUy), 

(18) 


where Ka and Ks are the absorption and scattering grey opaci¬ 
ties respectively and AnB = aT^ (a is the radiation constant) 
is the intensity of black body radiation for a gas with temper¬ 
ature Tg. We further define an effective radiation temperature 
Tr in the fluid frame as 

E = aT^. (10) 

The Comptonization component of Equation (8), under the 
“blackbody” approximation (see Sadowski & Narayan 2015) 
is. 


where di{ ^f^pu^) = 0 (eq. 14) due to the assumption of a 
stationary state. We move the metric determinant out of the 
derivatives and combine this with Equations (17) and (18) to 
obtain. 
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where is the electron mass. 

In the orthonormal fluid frame, equation (9) corresponds 
to the more intuitive 

G‘ = Kap{E-47:B), (12) 

G‘=iKa + Ks)pF\ (13) 

where F is the radiation flux. Equation (12) describes the 
energy transfer rate resulting from absorption and emission 
while Equation (13) describes the rate of change of the mo¬ 
mentum because of absorptions and scatterings. 

Equations (4)-(6) can be rewritten in a coordinate basis by 
separating the time dimension from the spatial coordinates 
yielding, 

dt(pu‘) + dii ^pgpu') = 0 (14) 

diiTi) + dii yFsK) = + yTgGy (15) 

d,iR[) + dii V=gO = ^FgR'‘,Yi - y^gGy (16) 

where is the metric determinant and are the 

Christoff el symbols. 


2.2. Decomposition 

To decompose the forces acting on the gas, we rewrite 
Equation (15) assuming a stationary state {dt = 0) and de¬ 
composing the stress-energy tensor Ty from Equation (7), 


Here, the left hand side corresponds to the convective deriva¬ 
tive of the gas, while the terms on the right hand side corre¬ 
spond to different forces acting on the gas. We define them 
as follows. 


/v,thermal ~ ^y(/^g) 

W 


/> 


y,magnetic 


--dy(\b^\+Uimy) 
w \2 w 
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f — 

7y,metric ~ ^ y/c 

W 

/y,radiation — Gy 

W 


-pu^Uy 
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/y,enthalpy — P)^ ) • 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


Equation (20) describes the thermal pressure force. The 
thermal force is dependent on the gradient of the gas pressure 
dy(pg). This indicates that the force is due to the internal 
pressure the gas exerts on itself (or rather, the pressure the 
gas as a whole exerts on individual particles). 

Equation (21) is composed of two terms, the first which 
shall be referred to as the magnetic pressure and the second 
as the magnetic tension. Like the thermal force, the magnetic 
pressure is a gradient force. This force arises from the density 
of the magnetic field lines. Magnetic tension, on the other 
hand, is the force exerted on particles as bent magnetic field 
lines straighten out. 

Equation (22) describes the forces resulting from the space 
time curvature. Most importantly, one of its components 
reflects the gravitational acceleration. Another component, 
reflects the (virtual from the point of view of a stationary 
observer) centrifugal force. The remaining terms account 
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mostly for the relativistic corrections due to space-time cur¬ 
vature and are negligible far from the BH. 

Equation (23) describes the force that radiation exerts on 
particles. Photons can either accelerate particles due to radi¬ 
ation pressure, or decelerated them because of the radiative 
drag, as discussed above. 

Equation (24) is a relativistic correction resulting from the 
difference between the total system enthalpy w and the rest 
mass density p. 

We can furthermore decompose Equation (22) into the 
components driving the gravitational attraction and the (az¬ 
imuthal) centrifugal force. 


where. 


/v,metric 

~ /v,gravity /v,centrifugal /v,residual, 

(25) 


jt 

f _ t T-'t CV 

7 v,gravity — A yt^r 

W 

(26) 


rji(p 

r _ ^ 

yv,centrifugal ~ ^ vd)’ 

W ^ 

(27) 


The remainder, /y,residual, together with Equation (24) reflects 
the relativistic correction terms, denoted by 


/v,correction — /v,residual /v,enthalpy- (^8) 

The forces can be expressed in an orthonormal basis (as¬ 
suming a diagonal metric of a non-rotating BH) through 

fv — fv! 

Eor a non-relativistic system where u^-\- p^-\- ^ p and 
V 1, Equations (20)-(28) reduce to the well known Newto¬ 
nian formulae. 


weak magnetic field lines. In the case of the radiative run, 
local thermal equilibrium between gas and radiation was as¬ 
sumed. The radiation was evolved adopting the Ml closure 
scheme. Eor both runs the BH spin was zero. 

The simulations assumed axisymmetry and ran in 2.5 di¬ 
mensions, that is, assuming axisymmetry but allowing non¬ 
zero azimuthal components. We prevented quick dissipation 
of the magnetic field in axisymmetrical simulations by ap¬ 
plying the mean-field dynamo of Sadowski et al. (2015). 

The simulations were run for an exceptional amount of 
time, reaching ~ 2,000 orbits at the innermost stable cir¬ 
cular orbit (~ 200, OOOGMbh/^^^)- This allowed for the out¬ 
flow to reach inflow/outflow equilibrium (defined as the re¬ 
gion where average properties do not chang with time) in 
relatively large domain, especially in the polar region, where 
the gas moves with large velocities. Applying the criterion 
from Narayan et al. (2012) we obtain in both simulations the 
equilibirium regions extending up to ~ 80rg at the equatorial 
plane and up to the computational box boundary at ~ lOOOrg 
at the axis. 

Eor this work we adopted the following forms of the ab¬ 
sorption and scattering opacities, 

Ka =6.4 • 10^^ cm^lg (36) 

Ks =0.34 cmVg, (37) 

where p and Tg are the gas density and temperature, respec¬ 
tively. Such a form of absorption opacities reflects only the 
bremsstrahlung absorption and emission, and, in particular, 
neglects line opacities. Therefore, the radiative simulation 
applies only to hot accretion flows near stellar mass BHs. 


//, thermal — 

(29) 

Amagnetic = + ^di{b'bi) 

(30) 

/r,gravity - 

(31) 

f - • fl 

Jr,centrifugal ~ Sm u . 

^ r sm ^ 

(32) 

centrifugal ~ ^ 6 , 

^ rsmO 

/irradiation ~ ~Gi 

P 

(33) 

(34) 

/i, correction — 0 - 

(35) 

3. NUMERICAL METHOD 


We analyzed two simulations performed with a general 
relativistic (GR) radiation magnetohydrodynamical (RMHD) 
code KORAL (Sadowski et al. 2013; Sadowski et al. 2014) 
- one purely MHD (optically thin, hdSOOaO) and the other 
including radiation field (optically thick, d300a0) and cor¬ 
responding to the accretion rate ~ lOMpdd- Both were ini¬ 
tialized as equilibrium torii threaded with multiple loops of 


4. RESULTS 
4.1. Flow properties 

The leftmost and center panels in Eigure 1 show the av¬ 
eraged (top set of panels) and instantaneous (bottom set) 
density distribution and velocity vectors of gas in the non- 
radiative^ and radiative simulation respectively. The gas is 
clearly concentrated near the equatorial plane but the disk 
thickness is large - the density scaleheight in both cases is 
close to HjR = 0.3. In case of the hydro simulation, the gas 
relatively smoothly fills the polar region. Eor the radiative 
simulation, on the other hand, the polar region has signifi¬ 
cantly lower density than the bulk of the disk. In that case, 
the density at the axis can be as much as 6 orders of mag¬ 
nitudes lower than at the equatorial plane. This fact results 
from radiation cleaning the funnel, as will be discussed be¬ 
low. 

The streamlines in the left and center panels show the av¬ 
erage velocity field. In the case of the hydro simulation the 

^ Purely magnetohydrodynamical (non-radiative) simulations are scale- 
free, i.e., they correspond to optically thin flow of any density. Therefore, 
the values of density given above the left panels are arbitrary. 
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Table 1. Properties of the analyzed simulations. 


Model 

BH spin 

Mbh 

Radiative 

Accretion rate 

hdSOOaO 


lOMo 

No 


dSOOaO 

= 0 

lOMo 

Yes 

~ lOMndd 


gas flows inward only deep inside the bulk of the disk. Far 
from the equatorial plane vertical motion is clear, leading to 
gas escaping from the disk. The velocity field in case of the 
radiative simulation is noticeably different. Very fast (much 
faster than in the hydro case) outflow is evident only in the 
polar region. At intermediate polar angles, the gas does not 
show clear outflowing pattern, what implies that the turbulent 
motion is strong and on average does not give strong outflow. 

The corresponding bottom panels show instantaneous gas 
density distribution and velocity field. The turbulent nature 
of the flow, resulting from the magnetorotational instability 
(Balbus & Hawley 1991), is evident in both cases. It is inter¬ 
esting to note that in the hydro simulation gas can temporarily 
fall on the BH, although on average it shows positive radial 
velocity. These episodes feature very low density gas which 
does not provide significant inward mass flux. 

The right panels in Figure 1 show the averaged (top) and 
instantaneous (bottom panel) radiation field in the radiative 
simulation. Colors denote magnitude of the radiative flux and 
streamlines show its direction. Deep inside the disk photons 
are trapped and advected on the BH (Sadowski & Narayan 
2015). At intermediate angles radiation diffuses out of the 
disk and contributes to the radiation in the funnel. The effi¬ 
ciency of accretion in radiation emitted into the optically thin 
funnel is not large - radiation carries only ~ 1 % of Mc^, what 
is ~ 5 times less than it would carry for a thin disk accret¬ 
ing on a non-rotating BH. The corresponding bottom panel 
shows instantaneous properties of the radiation field. Radia¬ 
tion escaping through the funnel is quite laminar, in contrast 
to turbulent disk interior, where photons are dragged with op¬ 
tically thick, turbulent gas, producing, on average, the pattern 
described above. 

4.2. Average forces 

We begin by calculating the average forces acting on gas 
at a given location. We averaged the simulations’ output on 
the go, i.e., accounting for states after every single step of 
time iteration. Each product appearing in Equations 20-24 
was averaged independently, e.g., 

r; - pu‘Uy \ {Tl) - {pU'Uy) 

-/---• ^Jo) 

W I (w) 

In such a way we are able to obtain the average forces act¬ 
ing at a given location, i.e., the forces that observer sitting at 
fixed coordinates would see acting on gas flowing by. There 
is no guarantee that every single parcel of gas will feel ex¬ 
actly these forces when following its own trajectory. How¬ 
ever, one can hope that forces obtained in this way will give 


qualitatively right picture. 

The results of this analysis for the optically thin (hdSOOaO) 
and thick (d300a0) simulations are shown in Figures 2 and 3, 
respectively. The background color denotes the logarithm of 
the magnetic to rest-mass energy densities ratio. This quan¬ 
tity measures how strong the magnetic field is, and in partic¬ 
ular, how much the relativistic enthalpy w = p + + + 

departs from the non-relativistic value w = p. 

Using Equations (20)-(24) we calculated radial and polar 
components of each force and transformed them to the or¬ 
thonormal frame. Acceleration induced by each force com¬ 
ponent are plotted with vectors. Blue arrows correspond to 
the gravitational force (Eq. 26), green to the thermal pres¬ 
sure gradient force (Eq. 20), purple to the magnetic forces 
(Eq. 21), red to the centrifugal force (Eq. 27), and cyan to 
the relativistic correction (Eq. 28). In the case of the opti¬ 
cally thick, radiative run (Eig. 3), an additional set of yellow 
arrows denote the radiative force (Eq. 23). In both plots the 
black arrows correspond to the sum of all the forces, i.e., to 
the total acceleration acting on gas. The larger is the length 
of this vector, the faster gas particles gain velocity. All the 
vectors are scaled in the same way, with the scaling factor 
chosen arbitrarily to fit the arrows in the panels. 

4.2.1. Optically thin disk 

Eigure 2 shows the average forces in the optically thin, 
purely magnetohydrodynamical simulation. The two forces 
which are most obvious are the gravitational pull (blue) and 
the centrifugal force (red arrows). Gravitational attraction 
always works towards the BH and its magnitude decreases 
(up to relativistic correction) with the square of radius. The 
centrifugal force acts only in the direction perpendicular to 
the axis of rotation. Its strength depends on the value of az¬ 
imuthal velocity and the distance from the axis. The largest 
centrifugal acceleration acts on particles close to the BH and 
located near the equatorial plane, inside the almost Keplerian 
disk. It is much less significant in case of slowly rotating gas 
in the polar region. 

Outside of the polar region gas moves in a turbulent way 
with relatively low velocities and the velocity changes it un¬ 
dergoes are much less violent then for a free-falling particle. 
It is therefore reasonable to expect that the net force acting on 
the gas in the bulk of disk would be much smaller in magni¬ 
tude than the gravitational and centrifugal forces. To satisfy 
this condition one has to provide another force, or forces, 
which will balance the two. In case of optically thin disk it 
is the thermal pressure force (green vectors) which acts this 
way. It provides the necessary vertical component which bal¬ 
ances the vertical component of gravitational pull. As a re¬ 
sult, the total acceleration (black arrows) is small in the bulk 
of the disk. 

Situation is different in the polar region where the gas is ac¬ 
celerated rapidly and reaches significant radial velocities. In 
this case, the magnitude of the net acceleration may be even 
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Figure 1. Leftmost and center panels show the logarithm of the gas density distribution (in CGS units, glcrn^) and the gas velocity field for 
the non-radiative and radiative simulations respectively. The width of the vectors are proportional to the velocity magnitude. The rightmost 
panels shows the logarithm of the magnitude of the radiative fiux (erglsicrn^) and the radiative flux vector field. The top panels correspond to 
averaged data. The bottom panels represent instantaneous snapshot data at an arbitrary chosen time in the simulation. 


comparable with the gravitational pull. In contrast to parti¬ 
cles inside the disk, the centrifugal force is no longer substan¬ 
tial. Instead, the magnetic force (purple arrows) affects sig¬ 
nificantly the acceleration. Together with the thermal force, 
it provides the vertical force that balances and overcomes the 
gravitational pull. For the gas in the innermost polar region, 
where the magnetic energy density becomes comparable with 
the rest-mass energy density, one should not neglect the rel¬ 
ativistic corrections (cyan arrows), which may reach magni¬ 
tudes comparable with the other forces. 

4.2.2. Optically thick disk 


Figure 3 shows average forces acting on gas in the optically 
thick simulation. The colors once again denote the strength 
of the magnetic field. Because of much lower density of gas 
in the polar region of this simulation than in the case of op¬ 
tically thin disk, the magnetic field now dominates over rest- 
mass energy and satisfies Ip > 1 in most of the funnel 
region. 

The properties of the gravitational pull (blue arrows) and 
the centrifugal force (red arrows) are quite similar to the de¬ 
scribed above - the former points always towards the BH, 
while the latter is perpendicular to the axis. The centrifugal 
force points outward everywhere but for the points located 
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Figure 2. Averaged data from the non-radiative simulation. The fig¬ 
ure shows the ratio \og{b^ jp) (shades of grey) and the vectorial force 
distribution at arbitrary points where the arrows indicate both direc¬ 
tion and magnitude of the forces. 
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Figure 3. Similar to Figure 2 but for the averaged data from the 
radiative simulation. 


in the funnel region. This surprising effect comes from the 
fact that the component of equation (27) changes sign 
in that region because of predominant magnetic field com¬ 
ponent b^b(j) (if we had defined this force as proportional to 
pu^U(p, instead of the force would always point outward, 


as in the non-relativistic limit). Similar effect was not present 
in the optically thin case because magnetic field in the polar 
region was much weaker than in the case described here. 

Relatively slow gas velocities in the bulk of the disk re¬ 
sults from weak acceleration. This is once again the reason 
why the forces deep in the disk balance each other and the 
resulting average net force (and acceleration) is very weak 
(in most cases the corresponding arrow is not visible in the 
plot). However, on the contrary to the optically thin case, it 
is not the thermal pressure gradient which balances the verti¬ 
cal component of the gravitational force. In case of optically 
thick, radiation-pressure dominated disk like the one simu¬ 
lated, it is the radiation pressure force (equation 23, yellow 
vectors in Figure 3) which provides necessary vertical force. 
In the optically thick disk it follows the gradient of radiative 
energy density and therefore points away of the equatorial 
plane towards the axis. In other words, the radiative force 
has replaced the thermal pressure force and is now support¬ 
ing the disk against gravity. 

The properties of the forces change significantly in the po¬ 
lar region. As shown in the left panels of Fig. 1, strong flux 
of radiation coming from the innermost region is streaming 
out of the system along the axis. Gas which enters this region 
is immediately swept up and with mildly-relativistic veloci¬ 
ties moves away from the BH through the funnel. This ra¬ 
diative acceleration is clearly seen in Fig. 3 - radiative force 
in the polar region has significantly larger magnitude than 
the gravitational force at a given location and always points 
away from the BH. However, the radiative force is not the 
only one pushing the gas away from the BH. The magnetic 
force (equation 21) and the relativistic correction (equation 
24) also point upward and reflect the magnetic acceleration 
of gas in the magnetic tower that develops when magnetic 
field lines brought into the innermost region are wound up by 
disk rotation and expand vertically in the polar region drag¬ 
ging the gas behind. This dragging effect, and the resulting 
magnetic acceleration, is more effective in the optically thick 
case than it was in the other simulation because of lower den¬ 
sity of the gas and larger magnetic to rest-mass energy ratio. 

4.3. Trajectory approach 

The averaged data that we based on in the previous sec¬ 
tion reflects average forces acting at a single point over the 
course of the simulation, and does not reflect forces affecting 
a particular moving particle that happens to pass by the point 
in question at the specific time. What is more, not all the 
gas crossing given location, that on average belongs to the 
outflow region, is guaranteed to escape to infinity. There is 
also gas which either moves (probably temporarily) inward 
or is not energetic enough (is bound) to escape the BH grav¬ 
itational pull. It is therefore not straightforward to extend 
the analysis done in the previous section, which was based 
on averaged disk properties, to every single gas parcel. The 
force balance may be significantly different for outflowing 
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and inflowing gas crossing the same location. 

We therefore traced a series of particle trajectories for each 
simulation. We chose the trajectories to represent the jet, 
wind and disk region of the outflow. We injected virtual 
test particles and tracked their trajectories by interpolating 
the velocity from snapshot data at time t, and thus obtaining 
the positions of particles at time t + Ar. This method holds 
for sufficiently small Ar intervals. We used snapshot data 
saved every Ar = IGM/c^ what gave accuracy good enough 
to track even the fastest gas moving with mildly-relativistic 
velocities v ~ 0.3c. For each trajectory, we then calculated 
the decomposed real-time forces along that trajectory, which 
then were used to determine the acceleration mechanisms of 
the outflow 

4.3.1. Optically thin disk 

The top panel in Figure 4 shows the four representative 
trajectories that we chose for the optically thin model. We 
will refer to these trajectories as HI to H4 (numbered left to 
right). First two (HI and H2) correspond to gas entering the 
funnel and escaping along the axis. The other two (H3 and 
H4) reflect typical trajectories of gas escaping in the wind 
region of the simulation. The colors along the trajectories 
reflect the velocity of the gas. As expected, gas in the funnel 
undergoes the strongest acceleration and reaches velocities of 
the order of 0.1c. The wind trajectories (H3 and H4) reflect 
much slower gas, moving on average with velocities ~ 0.01c. 

The four panels below the top one show the forces acting 
on gas along each of the trajectories. The arrows have the 
same meaning as in the previous plots. The numbers on the 
plots reflect the Bernoulli number of the gas (equation 3) at 
each location for which we calculated the forces. Positive 
values correspond to unbound gas. 

The jet trajectories (HI and H2) both accelerate and then 
decelerate while maintaining a growing Bernoulli number 
making them part of the real outflow. However, the main ini¬ 
tial acceleration mechanism of these two trajectories differ, 
as trajectory HI is initially accelerated by the magnetic force 
(magenta arrows) and trajectory H2 is initially accelerated by 
the thermal pressure force (green arrows). It is evident that 
for these trajectories the acceleration is not a continuous pro¬ 
cess - gas is rapidly accelerated in the innermost region, it 
enters the funnel with large velocity, but then the further ac¬ 
celeration ceases and the gravitational pull slows down the 
gas. There is no unique force which is responsible for this 
initial acceleration - it could be either magnetic or thermal 
force. 

Wind trajectories H3 and H4 both exhibit relatively low 
positive radial velocities and show negative but growing 
Bernoulli numbers. They are part of the turbulent outflow 
within our simulation boundaries, but the growing number 
indicates that they could be part of the real outflow. We do 
not observe a sudden acceleration, but rather a series of minor 
acceleration events giving the trajectories on average positive 


radial velocities. The main acceleration mechanisms for the 
trajectory closest to the funnel region (H3) and the one deeper 
in the disk (H4) are the magnetic and thermal forces, respec¬ 
tively, with the centrifugal force also providing significant 
component along the trajectory. 

4.3.2. Optically thick disk 

Figure 5 presents in a similar way four outflowing tra¬ 
jectories in the optically thick simulation. We will refer to 
them as trajectories R1 (leftmost) to R4 (rightmost). The top 
panel shows the trajectories on top of the average gas density 
distribution. The colors along them denote the gas veloc¬ 
ity. Trajectory R1 stands out as the one closest to the axis 
- the gas following this trajectory reaches the highest veloc¬ 
ity, V ~ 0.4c. The second leftmost trajectory, R2, follows 
roughly the boundary between the polar funnel and the disk. 
The corresponding gas is not accelerated so efficiently but 
still reaches velocity of the order of ~ 0.1c. The other two 
trajectories are deep inside the optically thick region and are 
much slower with gas not exceeding 0.02c. 

The four panels below the top one show forces acting along 
each of the trajectories. Gas which managed to enter the fun¬ 
nel region (reflected in path Rl) has very low density and is 
strongly accelerated by the radiation force (yellow arrows) 
which is responsible for effective and continuous accelera¬ 
tion. This acceleration results in a mildly-relativistic outflow 
of unbound gas along the axis. This outflow is characterized 
by lowest densities but may carry significant amount of ki¬ 
netic energy (S^dowski & Narayan 2015). 

The gas following the boundary of the funnel and the disk 
(trajectory R2) is also pushed forward predominantly by the 
radiation force. This time, however, the force is weaker and 
is pointing more towards the axis, not along it. This results 
from the fact that at this location the radiation diffuses out of 
optically thick disk following the gradient of radiative energy 
density. Radiative flux, together with the centrifugal force, is 
strong enough, however, to accelerate the gas and make it 
unbound. 

The wind trajectories (R3 and R4) correspond to gas flow¬ 
ing through surface layers of the disk. The radiative flux in 
the comoving frame of the gas, which results from photon 
diffusion, points towards the disk surface and so does the ra¬ 
diative force. The radiative and centrifugal forces are not 
strong enough, however, to efficiently accelerate the gas - it 
is just enough to move the gas slowly outward. The corre¬ 
sponding Bernoulli numbers are negative but growing and it 
is possible that ultimately the radiation will manage to make 
the gas unbound. The properties of the force balance along 
these two trajectories very well agree with the results of the 
analysis based on averaged data. 

5. DISCUSSION 

The outflow in the optically thin, non-radiative simulation 
is mainly driven by the magnetic and thermal forces. The 


Acceleration oe wind in BH accretion disks 


9 


100 



A 


Gravity 

Thermal 

Magnetic 

Centrifugal 

Rel. correction 

Total 


i-4 


--5 


-6 


1-7 










x/rg x/rg 


Figure 4. Trajectories for optically thin model hdSOOaO. The top 
panel shows two jet and two wind/disk particle trajectories (H1-H4, 
left to right), where the color is proportional to the speed of the 
particle in fractions of the speed of light plotted over the logarithm 
of gas density. The four bottom panels show the force decomposi¬ 
tion at arbitrarily chosen points for each individual trajectory. The 
length of the force vectors is proportional to the force magnitude, 
scaled by r^, and the width of the trajectory is proportional to the 
particle velocity (for exact value, see top panel). Bernoulli number. 
Be, along the trajectories is shown. 


Figures. Similar to Figure 4 but for the optically thick model 
dSOOaO and trajectories R1-R4 (left to right). The additional ra¬ 
diative force is denoted by yellow arrows. 

initial acceleration mechanism for the radiative simulation is 
the radiative force. In highly magnetized and the innermost 
regions in both models the relativistic correction (Eq. 28) is 
not negligible, and the standard Newtonian decomposition of 
forces would not be adequate. 

The jet trajectories (HI and H2, see Section 4.3.1) left the 
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simulation boundary decelerating (gravitational pull was pre¬ 
vailing) but with positive and still growing Bernoulli num¬ 
bers, making them part of the real outflow, while the wind 
trajectories (H3 and H4) of the non-radiative simulation were 
characterized by negative but still growing Bernoulli num¬ 
ber. The maximal velocities in the non-radiative model were 
approximately 0.10c. The radiative accretion flow revealed 
much larger velocities, approaching 0.40c for low-density 
gas in the jet region accelerated vertically by the uniform ra¬ 
diation flux. Both radiative jet trajectories we studied (R1 
and R2, see Section 4.3.2) left the simulation bounds with 
positive and growing Bernoulli parameters. As in the non- 
radiative model, the gas contained in the wind was bound in¬ 
side the domain, but the outward forces were still increasing 
its Bernoulli number. 

In case of the optically thin, non-radiative disks, the ini¬ 
tial acceleration seems to play the major part in generating 
real outflow, as the gas is accelerated rapidly in the inner re¬ 
gion and then propagates with forces near equilibrium. The 
trajectories in the jet region of the radiative simulation, on 
the contrary, are continuously exposed to radiation force and 
thus undergo continuous acceleration. 

The average forces acting on the gas at a given location 
reflect dominating forces acting on a gas parcel along its tra¬ 
jectory. However, the average forces do not allow to catch 
the exact timing and location of the short-lived but signifi¬ 
cant acceleration episodes. They also do not reflect the turbu¬ 
lent nature of the flow. The trajectory approach provides full 
information about the acceleration mechanism for traced tra¬ 
jectories, making it ideal for studying this phenomena. While 
it is not possible to trace every single gas trajectory, the av¬ 
eraged data serves as a reference of what one should expect 
given a large enough set of traced trajectories. 

We have compared an optically thin disk, corresponding 
to the radiatively inefficient regime and to the lowest accre¬ 
tion rates, with a model of optically thick disk accreting at 
a super-ciritcal rate of lOMEdd- To get full understanding 
of processes driving outflows in accretion disks, one would 
have to extend our work with a similar study covering wider 
range of accretion rates, including geometrically thin, sub- 
Eddington disks, and taking rotating BHs into account. Es¬ 
pecially in the case of spinning BHs one may expect that the 
extra energy extracted from the BH spin will significantly al¬ 
ter the acceleration mechanisms in the jet region. 

Similar modelling of forces acting on gas in accretion 
flows was done previously in Yuan et al. (2015) and Taka- 
hashi & Ohsuga (2015). The former work focused solely on 
optically thin disks, while the latter analyzed only a super¬ 
critical disk similar to our radiative model. Our conclusions 
are in general agreement with both these works. However, 
we calculated and decomposed the forces in fully relativistic 
way and directly compared driving mechanisms in optically 
thin and thick disks. 


6 . SUMMARY 

We have analyzed the force balance and the resulting ac¬ 
celeration of the outflow in two simulations of a BH accre¬ 
tion flow. Both were performed in general relativity with 
a state-of-the-art numerical methods and assumed zero BH 
spin and weak, non-saturated magnetic held. One model 
corresponded to radiatively inefficient, optically thin disk, 
while the other to a super-critical disk with accretion rate of 
lOMEdd- We found that: 

1. Gas acceleration: Driving gas out of the disk is not, in 
most cases, a continuous and simple process. There are 
usually only a couple of short-lived episodes which re¬ 
sult in significant outward acceleration of the gas. But 
for them, the gas velocity fluctuates, as does the net 
force acting on gas. An exception is the polar region 
in radiative disks, where gas is constantly pushed away 
by radiation pressure. 

2. Outflows in optically thin disks: Radiatively inefficient 
accretion produces outflows at large range of angles. 
The gas in the polar region is accelerated mostly by the 
magnetic forces and reaches velocities of the order of 
0.1c. At larger distance from the axis, it is the thermal 
pressure and centrifugal forces which compensate for 
gravity and push the gas out of the disk. The magnetic 
forces are not significant but for the most polar region. 
Therefore, the outflow outside of the polar region is not 
driven by the magnetocentrifugal acceleration. 

3. Outflows in optically thick disks: The fastest outflows 
in radiative disks take place near the axis, in the opti¬ 
cally thin funnel region. Gas can reach there mildly- 
relativistic velocities exceeding at times 0.4c. Such 
significant acceleration results mostly from the radi¬ 
ation pressure force pushing the gas along the axis. 
Magnetic acceleration is also important in the polar 
region. In the optically thick disk it is the radiation 
pressure, not the thermal pressure as in the other case, 
which balances gravity and pushes gas out of the disks, 
either into the funnel or radially outward. The mag¬ 
netic acceleration is once again not significant but for 
the most polar region. 

4. Relativistic corrections: We found that the relativistic 
correction that we have singled out from other forces 
(Eq. 24) contributes significantly in the regions close 
to the BH and where magnetic held energy is compara¬ 
ble or exceeds the rest-mass energy density. Therefore, 
force decomposition based on non-relativistic formu¬ 
lae would be inadequate in these regions. 

5. Average trajectory approach: Only the study of bal¬ 
ance of forces based on following particular gas tra¬ 
jectories is able to catch the non-uniform nature of gas 
acceleration. 
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